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Cl 

A  effective  contact  surface  between  two  phases. 
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p  pressure. 

in 
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q  heat  flux. 
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water  vapor. 
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IX 


1.  INTRODUCTION 

Metal  particle  ignition  and  burning  phenomena  have  been  used  for  several  purposes. 
In  solid  rockets  metal  ingredients  enhance  the  specific  impulse  and  stabilize  acoustic 
waves.  In  under-water  vehicles  metal  particle  burning  is  used  as  fuel  for  its  motor. 
Recently  explosive  enhancement  of  Under-Water  (U/W)  Shaped  Charge  (S/C)  by  means  of 
aluminum  particle  burning  in  water  has  been  examined,  (Ref.  1). 

In  the  case  of  stand-off  explosion  of  an  U/W  S/C,  a  diverging  shock  wave  propagates 
into  the  surrounding  water.  A  hypersonic  metal  particle  jet  is  moving  in  a  radial  direction, 
from  the  exDlosive  center  to  the  target.  As  a  result  of  severe  heat  transfer  conditions  at 
the  stagnation  point  ofthe  S/C  metal  jet,  some  particles  may  ignite.  Ignition  will  depend  on 
initial  particle  size,  environmental  condition  and  time  exposed  to  heating  fluxes.  Combus- 
tion will  enhance  the  pressure  driving  the  shock  wave  and  tear  the  target  vigorously.  Vari- 
ous combustion  parameters  affect  the  impulse  on  the  target  in  an  U/W  S/C.  Some  of  them 
are  :  heat  of  combustion  of  metal  with  steam,  initial  size  of  burning  particles,  flow-field  etc.. 
A  close  investigation  of  different  type  of  metals  burning  in  steam  was  carried  out  by 
Exploding  Metal  Wire  (EMW)  tests,  (Refs.  2-4  ).  The  purpose  of  these  tests  was  to  find  the 
best  available  metal,  which  due  to  its  burning,  will  improve  the  U/W  S/C  effectiveness. 
Although,  as  shown  later,  there  is  no  similarity  in  combustion  behavior  between  the  U/W 
S/C  and  EMW,  one  can  assume  that  the  most  effective  and  energetic  metal  found  in  the 
EMW  tests,  will  also  be  the  best  to  use  in  U/W   S/C. 

Ref.  1  describes  some  experiments  with  aluminum  shaped  charge  and  a  calculation  of 
impulse  improving  by  factor  of  1.2  due  to  complete  droplet  combustion.  Strong  transient 
effects  predominate  the  S/C  particles  heating  process;  therefore,  ignition  and  combustion 
must  be  carefully  examined.  In  some  cases  large  particles  would  not  ignite,  (Ref.  5),  and  in 
others  existing  burning  particles  may  be  extinguished.  The  problem  of  predicting  the 
amount  of  metal  particle,  varied  by  size,  that  contributes  by  combustion  to  the  explosive 
enhancement,  needs  to  be  solved  by  constructing  an  appropriate  model.  Before  postulating 
a  new  model,  we  need  to  understand  the  various  parameters  and  stages  in  the  general 
combustion  process.  We  will  discuss  the  metal  ignition  problem  of  U/W  S/C  in  water  vapor 
environment.  Then  we  will  examine  the  the  governing  parameters  in  the  fast  combustion 
process  and  construct  an  appropriate  model  by  starting  from  the  differential  equations  that 
describe  the  burning  process. 


2.  GENERAL  STAGES  IN  COMBUSTION 

There  are  three  distinct  time  intervals  in  the  combustion  process:  ignition, 
burning  and  cooling.  The  ignition  stage  includes  heating  of  the  metal  surface  to  a  certain 
level  of  temperature  where  slow  chemical  reaction  will  be  replaced  by  near  infinite  reac- 
tion rates.  As  a  result  of  rising  wall  temperature,  surface  reaction  rate  increases  and  thus 
the  oxide  layer  grows.  Heat  is  generated  by  the  reaction,  and  the  progress  of  the  burning 
process  depends  on  positive  supply  of  heat  to  the  metal.  If  heat  loss  rate  is  larger  than  the 
rate  of  heat  released  by  combustion,  the  burning  will  be  extinguished.  The  equation  that 
describes  the  heat  transfer  balance  is  given  in  Ref.  5  as  the  following: 


(1./3.)pMc      TdlldX  =  PoxHfuaS/fl  -  (k/2r)Nu(T-T)  (2-1) 

Ms 


Here  5  is  the  thickness  of  the  oxide  layer  and  can  be  presented  by  an  Arrhenius  type 
equation.  Eq.  2-2  from  Ref.  5  describes  the  kinetic  law  of  oxidation. 


db/tt  =  (B  coxm/5n)  Exp(-E/(RT))  (2-2) 


Spontaneous  high  heating  rate,  like  an  explosion  or  fast  electrical  heating,  may  cause  igni- 
tion with  a  very  small  growth  of  oxide  layer.  Fig.  1  illustrates  the  difference  in  temperature 
rise  between  the  fast  and  slow  heating  process.  In  the  spontaneous  case,  fast  heating  rate 
raises  temperature  high  above  the  ignition  point  and  then,  when  the  combustion  is  starting, 
the  temperature  drops  below  the  level  of  adiabatic  flame  temperature.  In  contrast,  com- 
bustion due  to  low  heating  will  occur  when  the  metal  surface  temperature  reaches  the 
ignition  temperature. 
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Fig.  1  Temperature  rise  in  fast  and  low  heating  rates 
Some  materials  such  as  aluminum  and  magnesium  have  a  Metal  Boiling  Temperature 
(MBT)  lower  than  its  Oxide  Boiling  Temperature  (OBT).  In  such  volatile  material,  combus- 
tion begins  just  after  the  metal  surface  is  exposed  to  the  surrounding  atmosphere.  Expo- 
sure of  metal  surface  occur  due  to  surface-tension-induced  movement  of  liquid  oxide  layer 
gliding  over  the  liquid  metal  or  by  blow-off  of  the  oxide  layer  as  happens  in  the  exploding 
wire  case.  Vapor  phase  reaction  will  take  place  in  the  flame  zone  as  long  as  heat  losses 
are  lower  than  heat  generated  by  the  combustion  process.  Glassman  stated  in  Ref.  6,  that 
the  maximum  temperature  in  the  combustion  zone  is  bounded  by  the  oxide  evaporation 
temperature. 


Non-volatile  metals  such  zirconium  or  tantalum  has  MBT  greater  than  OBT,  and  since 
the  combustion  temperature  is  lower  than  OBT,  the  metal  would  not  be  able  to  evaporate; 
therefore,  combustion  will  occur  on  the  metal  surface,  and  the  oxide  product,  near  the 
reaction  zone,  will  appear  on  the  gas  phase  only.  Table  1,  from  Ref.  7,  presents  some  data 
on  MBT  and  OBT  for  several  useful  metals. 


Metal 

Melting 
Pt.,°C 

Boiling 
Pt.,°C 

Oxide 

Melting 
Pt.,°C 

Boiling 
Pt,,°C 

Adiabatic 
Flame  Temp. ,°C 
(Equilibrium) 

Na  (Sodium) 

.97.8 

883 

Na20 

— 

1770   - 

i827 

Mg  (Magnesium) 

650 

1103 

MgO 

2800 

3100 

3077 

Ca  (Calcium) 

850 

1690 

CaO 

2590 

3500 

3527 

Al  (Aluminum) 

660 

2500 

A12°3 

2050 

3500 

3527 

Ti  (Titanium) 

1730 

3260 

Ti203 

2100 

3000 

3027 

Zr  (Zirconium) 

1845 

4900 

ZrO 

2690 

4300 

4527 

TABLE  1  Physical  properties  of  some  metals  and  their  oxides  (Ref.  7) 


The  composition  of  the  environment  has  a  strong  influence  on  the  combustion  process. 
Static  pressure  as  well  as  oxidizer  concentration  and  type  of  ingredients  determine  the  rate 
of  combustion.    Burning  of  an  aluminum  particle  in  free  air  is  governed  by  vapor  phase 


reactions,  whereas  in  rocket  motors  various  species  such  as  CO2  add  some  surface  reac- 
tions to  the  vapor  phase  burning  droplet  (Ref.  8).  Volatile  metal  burning  is  controlled  mainly 
by  vapor  phase  reactions  and  produce  large  quantities  of  small  condensed  oxide  particles, 
often  called  smoke  oxide  cloud.  The  small  condensed  particles  (less  than  2  micron,  Ref.  9) 
are  produced  in  the  flame  zone  and  convected  away  by  the  main  flow.  If  the  flame  zone  is 
close  enough  to  its  moving  burning  droplet,  smoke  oxide  particle  may  impinge  on  the 
metal  surface,  deposit  on  it  and  change  its  mass  and  momentum.  For  large  particles  the 
effect  of  smoke  oxide  radiation  may  not  be  neglected.  King,  in  Ref.  8  neglected  radiation 
effects  and  therefore  limited  his  solution  to  a  maximum  of  30  micron  aluminum  particle 
size,  and  found  good  results  compared  with  experimental  tests  in  wet  oxidizer,  Ref.  16. 
Wood,  Ref.  10  observed  strong  radiation  effect  on  burning  magnesium  particles  greater 
than  80  micron.  Heat  flux  radiation,  emitted  from  the  condensed  cloud,  cools  the  reaction 
zone  and  heats  the  evaporating  metal  surface.  A  schematic  description  of  the  various 
zones  in  a  burning  moving  aluminum  particle,  is  given  in  Fig.  2,  taken  from  Ref.  9. 
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Fig.  2  General  nature  of  the  combustion  zone   of  a  volatile 
metal  droplet  burning  in  a  mild  convective  flow  (  Ref.  9) 


3.  ESTIMATING  THE  MAXIMUM  PARTICLE  SIZE  FOR  IGNITION  IN  AN  U/W  S/C 

The  first  step  for  estimating  impulse  enhancement  of  an  U/W  S/C  by  means  of  alumi- 
num combustion  is  to  assure  that  the  moving  particles  in  the  S/C  jet  do  ignite  under  the 
special  conditions. 

3.1  Conditions  and  Method  of  Estimation 

In  order  to  make  a  quick  estimation  whether  an  aluminum  particle  can  reach  the  burn- 
ing phase  condition,  while  moving  at  high  speed  inside  a  shaped  charge  jet  of  particles 
cloud  at  an  underwater  explosion  conditions,  we  evaluate  the  time  needed  by  a  single  par- 
ticle to  increase  its  temperature  to  the  melting  point  and  become  completely  molten.  If 
that  time,  assuming  severe  conditions,  is  larger  than  the  transit  time  from  the  stand-off 
explosion  to  the  target  then  it  is  impossible  for  burning  conditions. 

Time  dependent  conditions  are  not  taken   into  account,  and  a  simple  calculation   is 
being  done  by  assuming  an  equilibrium  state  between  water  and  steam,  which  means  that 
the  cavity  is  full  with  saturated  steam.    These  conditions  are    the  most  severe  boundary 
condition  the  shaped  charge  particles  can  meet. 
Other  assumptions  made  here,  are: 

-The  surrounding  pressure  is  equal  to  the  hydrostatic  water   pressure. 

-All  particles  meet  the  same  conditions  and  donot  exchange   heat  betwe  en  themselves. 

-The  particle  spherical  shape  does  not  change. 

-Particle  velocity  remains  constant  and  equal  to  7000  m/sec. 

-Internal  temperature  distribution  is  ignored. 

3.2  Equations 

For  the  solid  state  particle  the  equation  of  lump  heating  is: 
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>.  (T  -T)-g(e  T    -  a  T    )  =  -pM  c       (dT/dt)  (3-1 ; 

cx    r      '      x  w  9    9  Ms    p.,   x  ' 


Here    X,  ,T  ,T    and  a    will  be  calculated  according  to  the  steam  conditions  and  the  particle 

c'    g'    r  g  ** 

velocity.   The  second  equation  shows  the  change  in  particle  radius  during  melting  time 


Here    X.  ,T  ,T    and  a    will  be  calculated  according  to  the  steam  conditions  and  the  particle 

c     g     r  g  ~  ~ 

velocity.   The  second  equation  shows  the  change  in  particle  radius  during  melting  time 


VTr  •  TnJ-a(EwTMm     "  agTg  )    "    "^Ms  "  ^iKsip^  &j ft  (3-2) 

3.3  Calculation  of  the  Coefficients 

a.  Saturated  Steam  Condition  at  3  atm. 

T      =  407  K   ;     p       =1.68  kgm/m3 

wv  rwv  ° 

Rw=  848/18=  47.11  kgfm/kgm*  K 
c      /  R     =  4.12  c 

p  wv  p 

wv  wv 

v  =  (1./(1  -  c      /  R      )  =  1.32 

'  x         v  p  WV   ' 

wv 


The  speed  of  sound,  a: 

a=  VryRTgo  =  /y1.32*47.1 1*9.81*407  =  498.38  m/sec. 

Viscosity  at  T  :     p  =  15.*10"      kgm/m*sec 

Thermal  conductivity  at  T  :    k  =  6.5*10"6  kcal/m*sec*  K 

Radiation  absorbance,  assuming  equivalent  length  of  2  meter  and  a  pressure  of  3  atm: 

«g(3atm-2J=    °-7 
b.  Recovery  Temperature  and  Stagnation  Conditions. 


M  =  u/a  =  (7000.  /  498.38  )  =  14. 
Tr  =  T(1+  (Y-1)rM2)/2) 
where  r  is  the  recovery  factor.  Assume  r=  0.89,  for  turbulent  flow. 
Tr  =407.(1+  0.16*0.89*14*)=  11766  K 


This  is  a  very  high  temperature,  and  we  can,  therefore,  expect  dissociation  and  an 
improvement  in  the  thermal  conductivity  and  viscosity.  The  density,  viscosity  and  Prandtl 
number  at  the  stagnation  condition  are: 

p°  =(  p/  (R^TJ  =  3.  104/(47.iri176.)  =  0.054  kgm/m3 
/  =40/10"6   kgm/m*sec 
Pr°  =  0.88 

c.  Aluminum  Properties. 
pMs  =  2700  kgm/m3 

pMf  =  2380  kgm/m 

Pms  "  PMf  =  320  k9m/m 
TM    =  932  K 

hMsf  =  94.45  kcal/Kgm 

c        =  0.23  kcal/kgm*  k 

PMs 

ks    =  0.05  kcal/rrTsec*  k 
*w    -  °-2 

d.  Stagnation  Heat  Transfer  Coefficient. 

The  heat  transfer  coefficient  is  defined  as: 

X.    =     0B   Ue   C        St 
Pe 

St    is  the  Stanton  number  for  stagnation  point  flow  around  a  sphere 

ofi  ofi         o    a6     o       as       o       o^ 

St  =  0.763  (pD/u  )      (ReD)       (Pr  )      (p  /p)      (p  /p)  P 

For  Mach  number  greater  than  10,  DJVu  =  1. 

In  order  estimate  the  shortest  time  for  melting,  we  will  assume  a  small  particle  diameter 
of  the  size  of  20  micron,  the  appropriate  Rep  number  will  then  be: 
ReD  =1.68*7000*20*10"6/15*10~6    =15680. 

The  ratio  of  density  at  stagnation  and  free  stream  and  the  ratio  of  viscosities  at  these  two 


conditions,  are 
p  /p  =0.054, 
The  heat  transfer  coefficient  is; 


p°/p  =0.054/1.68=  .032     ;  p°/p  =40/15=  2.67 


X  =  1.68*7000*.453*0.763*(15680)"(q,5)*0.88"a6#(0.32*2.67)O^=  10.23  kcal/m'sec  K 


c 


3.4  Comparison  Between  Forced  Convection  and  Radiation  Heating 

The  convection  term  is  defined  as: 

qc  =  Xc(Tr-T) 
Assume  particle  initial  temperature  of  500  K 

qc  =  10.23*(  11766-  500.)=  1.153*10*    kcal/m2*sec 
Maximum  radiation  heat  losses  that  can  occur  at  the  metal  surface  are  at 
the  highest  temperature,  T  =  Tm 

qr  =<y(£wTMn|4-o gT  4)  =  1.36*10""(0.2*9324    .  -  0.7*4074)  =1.46  kcal/m'  'sec 

As  qc  is  much  greater  than  qr   we  may  neglect  radiation  heat  losses. 

3.5  Diameter  Limit  for  Ignoring  Temperature  Distribution 

Temperature  gradient  inside  matter  can  be  neglected  when  k  /X'r  is  greater  than  6.  In 
that  case  :  0.05/10. 23*r  is  greater  than  6.  The  condition  is  valid  for  every  particle  diameter 
less  than  1630.  micron. 

3.6  Calculation  of  Heating  Time 

The  calculation  of  the  heating  time  needs  the  solution  of  the  two  differential  equations 
(3-1)  and  (3-2)  dropping  the  radiation  term.  Without  the  radiation  form  the  solution  is  easy  at 
two  different  periods  of  time. 

a.  Conduction  in  the  Solid  State  Time  Interval. 

Equation  (3-1)  can  be  written  for  a  non-dimensional  variable   0,  as  follows: 

(d0/dt)  -  t"1©  =  0  (3-3) 

Where:  0  =  (TMm-  T.)/(  T  -T.)  ,T.    is  the  initial  temperature,  and  t  is  the  time  constant  and  it 
is  equal  to  :    pu  c      r    ,  3*X 

Ms    p„     po/  c 

MS 

The  solution  for  equation  (3-3)  is: 
0=1-  exp(-t/t) 

If  tm  is  defined  as  the  heating  time  for  a  particle  to  reach  the  nondimensional  melting 
temperature  0  ,  we  can  find,  from  the  solution  of  equation  (3-3),  the  ratio  of  tm/r  as  a 
function  of  the  volumeteric  heat  capacity  pu  c       and  the  heat  transfer  coefficient  "k 

Ms 
Vpo   =   -(PmsCpMs/<3)c>  I  M  1-0„>]  <3"4) 


b.  The  Melting  Time. 

Forming   the   integration    Eq.(3-2),    by   neglecting   the   radiation   term   and    introducing   the 
boundary  condition  of  %  =  1  when  t  =  t  ,  yields  : 

MTr  "TmJ  (*  "  tj    -   -  (hMsfrpo  /  3)(PMs  -  PMf)  &     '  1>  (3"5> 

End  of  melting  occurs  when  0  =  0  and  at  that  time  t  =  tf 


L  *  t„  +  [hM  j      (  pM  -pM.)  /  (  31  (  T    -  TM  )]  (3-6) 

t  m  L    Msf  po      *ns  rMf'       v        cv      r         Mm/J  v        ; 


Using  equation  (3-4)  for  t  ,  we  will  find  a  general  expression  for  the  minimum  complete 
melting  time  tf 


t,/r      =  (1/3X  )  [  -  pM  c      £n(  1  -  0  )  +  (  pM   -  pM.)  hM  ,  /  (  T    -  TM  )  ]  (3-7) 

f     po  v  c'  l      rMs   p„       v  nr  v  KMs      r\W'      Msf      v      r         Mm'  J  v        ' 

Ms 


The  initial  temperature  T.   of  the  particle,  at  its  initial  trajectory  point,  is  determined  by 
conservation  of  energy.   We  will  assume  that  T.  =  500  K,  and  t^r     will  be  found  by: 

0m=  (932. -500.)  /  (11766.-500.)    =    0.0383 

and  inserting  aluminum  metal  properties  into    Eq.  (3-7),  yields: 


t.  /  r      =  0.0874  -  20.23*£n(1-©  )  (3-8) 

f        po  m' 
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By  introducing  ©m=  0.0383,  into  Eq.  (3-8)  we  determine  tf/r  equal  66.2  sec/m  .  Let  us 
assume  a  traveling  time  of  200  microseconds  for  a  particle  to  reach  the  target.  Now  we  can 
find,  by  substituting  tf  =200  microseconds,  that  just  6  micron  particles  or  smaller  will 
completely  melt  and  particles  with  larger  diameter  will  have  no  chance  to  reach  the  igni- 
tion point. 
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4.  EXPERIMENTAL  RESULTS  IN  WATER  VAPOR  ENVIRONMENT 

As  mentioned  in  chapter  2,  combustion  behaviour  varies  with  the  change  of  the  envi- 
ronmental composition.  Generally,  the  burning  rate  depends  on  the  oxidizer  concentration 
(  Refs.  7,16  ),  therefore  combustion  in  water  vapor  is  faster  than  in  air  at  the  same  static 
environmental  pressure  and  temperature.  On  the  other  hand,  ignition  in  water  vapor  is 
more  difficult  than  in  air.  The  oxide  layer  deposit  in  wet  environments  is  more  protective 
than  in  dry  environments  and  therefore  ignition  temperature  in  water  vapor  is  1700  K  (  Ref. 
8),  whereas  in  dry  air,  about  1000  K.  The  predominant  reaction  of  aluminum  in  water  has 
the  form  of: 


2  Al  +  3  H20  =  AI2O3  +  3  H2 


Wilson  et  al.  (  Ref.  11  )  found  the  reaction  rate  constants  for  the  above  reaction. 


K  =  4.916  10S  Exp[  -73.5/(1. 99T)]         (4-1) 


Ref.  11  also  mentioned  that  the  condensed  oxide  smoke  is  carried  away  faster  in  moist 
steam  than  in  dry  air.  Ref.  7  found  that  hydrogen  attacks  the  oxide  deposition  layer  causing 
porosity  and,  therefore,  greater  vapor  metal  fuel  consumption. 

Experiments  with  exploding  wire  in  steam  were  carried  out  by  Kol  et  al.,  Refs.  2-4.  Fig. 
3,  from  Ref.  4,  shows  an  example  of  magnesium  exploding  wire  and  its  different  main  com- 
bustion regions.  The  white  brightest  sphere  at  the  center  is  the  plasma  core.  Its  tempera- 
tures high  above  6000  K,  and  it  contains  partially  ionized  metal  vapor.  The  outer  spherical 
shell  consists  of  metal  vapor,  where  a  temperature  higher  than  the  OBT  is  maintained.  On 
the  outer  surface,  the  metal  vapors  meet  the  surrounding  steam,  and  vapor  phase  reac- 
tions occur. 


Since  radiative    and  conductive  heat  fluxes  from  the  plasma  core  to  the  metal  vapor  shell 
are  high,  metal  vapor  temperature  is  expected  to  be  much  higher  than  the  OBT.  Because  of 
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electrode 
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Plasna     -     white 

nagnesiun  vapor  -greei 
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Fig.  3  Magnesium  exploding  wire  -  region  of  combustion  in  steam 
(from  Ref.  4) 
these  enormous  heat  fluxes  and  high  temperatures,  the  system  can  not  be  considered  as 
an  adiabatic  system,  and,  hence,  local  oxide  product  may  still  appear  on  the  gas  phase. 
Results  of  measurements  presented  in  Ref.  4  show  that  up  to  a  Normalized  Distance  (ND) 
of  120  from  the  plasma  core  center,  the  cloud  and  its  streaks  are  blue  and  the  temperature 
is  above  the  OBT.  Beyond  this  distance,  color  is  changing  and  temperatures,  detected  by 
a  two  color  photo-pyrometry  technique,  indicate  condensation  of  the  oxide  gas  species. 
Overall  maximum  resolution  of  the  measurement  device  is  2  micron,  and,  therefore,  it 
could  not  resolve  the  small  condensed  oxide  particles  in  the  streak.  As  a  consequence  of 
the  above  measurement  results,  the  author  concludes  that  the  cloud  around  the  exploding 
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wire  consists  of  plasma,  gas  phase  reactant  and  gas  products.  Due  to  high  temperatures 
in  that  region,  metal  vapor  and  oxide  product  condensation  is  not  possible  in  the  corona 
sphere  close  to  the  exploding  center.  Oxide  product  condensation  may  take  place  just 
after  expanding  away  from  the  explosion  center.  As  a  result  of  the  wire  explosion,  some  of 
its  molten  parts  are  blowing-off  as  metal  droplets  and  appear  beyond  the  corona's  streaks. 
Fig.  4,  which  presents  a  schematic  model  for  emission  from  the  metal  droplet,  is  based  on 
a  color  slide  picture  taken  in  a  magnesium  wire  explosion  test.  In  some  cases,  one  can 
assume  that  hot  droplets  emitted  from  the  hottest  melted  wire  may  ignite,  but  most  of  them 
do  not.  They  simply  cool  and  condense  along  their  trajectory  in  steam.  Ref.  2  describes 

2 

measurement  results  of  the  high  particle  deceleration  (of  the  order  of  100000.  m/s  ), 
caused  by  growth  of  liquid  alumina  on  the  surface  of  moving  burning  metal  droplets.  In 
most  deceleration  measurements,  in  Ref.  2,  the  solid  oxide  particles  were  the  product  of 
abrupt  solidification  of  molten  super  cooled  droplets. 

The  various  cases,  described  above,  present  the  complexity  of  the  phenomena  of 
metal  combustion  under  various  environmental  and  initial  conditions.  In  order  to  have  bet- 
ter estimates,  there  is  a  need  for  a  calculating  tool.  The  next  chapter  is  devoted  to 
describing  the  transient  combustion. 
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Fig.  4  Regions  surrounding  a  metal  exploding  wire 


5.  DISCUSSION  OF  PARAMETERS  FOR  MODELING  THE  METAL  COMBUSTION 

Exploding  wire  and  underwater  shaped  charge  metal  jet  are  just  two  instances  of  a 
very  fast,  non-steady  phenomena.  Even  metal  particles  ignition  and  combustion  in  a  rocket 
motor  environment  cannot  always  be  considered  as  a  quasi-steady  process.  A  precise 
prediction  of  ignition  time,  temperatures  and  species  existing  in  the  various  combustion 
regions  can    be  done  only  if  transient  effects  are  included. 
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All  quasi-steady  combustion  models  (Refs.  6,8,12  )  assume  Lewis  Number  unity,  which 
physically  means  that  the  evaporating  metal  mass  flux  is  independent  of  the  gas  phase 
reaction  rate.  That  assumption  may  not  be  valid  in  any  of  the  following  circumstances: 
Reactions  are  kinetically  controlled  ,  strong  mass  convection  effect  and  dissociation  caused 
by  high  temperature.  Those  chemical  and  flow  conditions  may  certainly  be  present  at  the 
beginning  and  the  end  of  the  combustion  process  of  a  moving  metal  droplet. 

The  fact  that  the  moving  burning  metal  droplet  accumulates  part  of  its  oxide  product 
induces  a  dramatic  changes  in  its  velocity,  Ref.  2.  A  physical  model  will  be  required  to 
describe  the  evolution  in  time  of  the  various  parameters.  In  order  to  achieve  that  model, 
solution  of  the  flow  field,  mass  and  heat  transfer  equations  is  needed.  It  is  also  important 
to  estimate  the  radiation  and  the  effect  of  deposition  of  smoke  oxide  cloud  particles  on  the 
burning  metal  droplet.  It  is  necessary  to  account  for  the  flow  as  two-phase  flow  together 
with  radiation  terms.  Radiation  fluxes  from  the  condensed  oxide  causes  heat  losses  from 
the  reaction  zone  and  heats  the  burning  surface.  As  a  consequence  of  two-phase  flow  heat 
and  momentum  are  delivered  from  one  phase  to  the  other,  and  convection  terms  between 
two  phases  change  temperatures  and  velocities.  Therefore,  these  effects  should  not  be 
ignored. 

For  the  same  kind  of  metal  and  environmental  conditions,  differences  in  results 
between  the  various  1-D,  quasi-steady  models  (  Refs.  6,8,12  )  and  unsteady,  2-D  two-phase 
flow  model  are  expected.  The  validity  of  including  the  volatile  and  non-volatile  metals  in 
one  model  has  not  yet  been  demonstrated.  The  idea  of  joining  them  into  one  model  is 
based  under  the  following  thermodynamic  facts.  The  difference  between  these  two  types 
of  metal  is  mainly  in  the  reaction  zone  location  and  behavior  of  oxide  product.  By  solving 
the  coupled  conservation  equations,  one  can  find  the  distribution  of  enthalpy  and  vapor 
concentration  distribution  in  time  and  space.  But  metal  surface  temperature  and  location 
are  yet  not  known.  By  writing  the  mixture  total  enthalpy  expression  and  including  terms  of 
species  concentration,  heat  of  formation,  metal  phase-change  enthalpy,  and  by  employing 
the  parameter  of  vapor  metal  quality  and  the  thermodynamic  relationship  of  wall  tempera- 
ture and  metal-vapor  concentration,  one  can  eliminate  the  instantaneous  surface  location 
and  its  temperature.  In  the  case  of  non-volatile  metal  combustion,  the  thermodynamic 
dependence  of  surface  metal  temperature  on  vapor  concentration  along  with  solution  of 
conservation  equations  will  result  in  the  highest  temperature  being  on  the  burning  metal 
surface.  Whereas  in  the  case  of  volatile  metal  burning,  the  point  of  highest  temperature 
may  be  anywhere  in  the  field.  Oxide  condensation  takes  place  whenever  the  local  gas 
temperature  is  lower  than  OBT.  Since  MBT  of  a  non-volatile  metal  is  greater  than  OBT, 
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oxide  condensation  may  appear  at  a  large  stand  off  from  the  metal  surface.  In  contrast,  for 
volatile  metal,  condensation  of  oxide  may  be  on  the  metal  surface.  The  above  discussion 
shows  that  by  applying  the  thermodynamic  phenomena  description  and  the  conservation 
equation,  a  general  transient  model  can  be  formulated  for  volatile  as  well  as  for  non-vola- 
tile metal  combustion.  The  next  step  will  describe  the  governing  equations  of  such  a 
model. 
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6.  GOVERNING  EQUATIONS 

The  various  cases  of  exploding  wire  or  droplet  metal  combustion  in  steam  can  be 
solved  by  one  general  model,  which  includes  the  transient,  two-dimensional,  asymmetrical 
effects.  For  the  sake  of  simplicity,  since  we  are  interested  in  investigating  the  transient 
two  phase  flow  and  radiation  effects,  we  will  first  consider  the  problem  in  transient  one-di- 
mensional, spherical  configuration. 

6.1  Energy  Equation  for  Gas  and  Liquid  Metal 

The  equation  will  be  written  for  enthalpy  and  will  include  in  one  region  (subscript  1) 
the  various  gas  species  and  liquid  metal  as  well.  Properties  at  each  point  will  be  pre- 
sented as  dependent  on  species  mass  fraction  and  local  temperature.  Assuming  constant 
static  pressure  in  the  entire  region,  the  equation  in  the  radial  direction  will  be  given  as  fol- 
lows: 


d(piht)/dt  =  r~2d[r2Ki(5T,/dr)]/dt  +  q"i  -  Pluid(h,/dr)  (6-1) 


Where  the  mixture  enthalpy  is  given  by: 


h,  =  Lh.Ci  (6-2) 


C:  is  the  mass  fraction  of  species  i,  and  h:  is  given  by: 


h.  =  h  .  +  jc    dT  @  i  =  Metal, H20,H2  (6-3) 


Metal    enthalpy,    hM    ,    will    account   for    quality    of   metal    vapor,    x,    and    phase    change 
enthalpy,hMfg 


hM  -  hMg  +  (1  "  x)hMf9  M 


The  heat  source,  Q  1     is  the  sum  of  other  various  sources,  Q  1 


a.  Heat  released  by  metal  combustion 


i 


Qn  -  wMH,..  (6-5) 


M     fu 


b.  Heat  released  by  condensation  and/or  solidification  of  oxide  product  : 


Qi2  =  wh,         @T=T.  (6-6a) 

ox    oxfg         v-'  oxb  v  ' 


Q  13  =  w      h__       @  T  =  T    .  (6-6b) 


ox    oxfs  ^-^  oxm 


Two  assumptions  are  made  here.  The  first  is  based  on  the  idea  that  oxide  phase 
change  will  occur  if  local  temperature  is  lower  then  the  OBT  or  lower  then  Oxide  Melting 
Temperature  (OMT).  The  other  assumption  is  based  on  the  fact  that  mass  transfer  between 
phases  is  a  comparatively  fast  process. 
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c.  Heat  transfer  from  one  phase  to  the  other. 


Q  15  =  (  h  A  /V  )(Tt-T2)  (6-7) 


d.  Radiation  heat  source  due  to  cloud  particles, 
q       and  q  rw  are  the  heat  radiation  sources,  located  at  the  cloud  and  at  the  wall  respec- 

c  „ 

tively.    Later-on,  we  will  develop  those  two  radiation  terms.    In  general,  q      stands  for  radi- 
ation flux,  and  and  the  radiation  source  will  be: 


Q  M  =  q  r/l  (6-8) 


The  solution  of  the  energy  equation  together  with  the  conservation  equations  will 
produce  an  enthalpy  distribution  for  each  time  step.  Temperature  and  metal  metal  vapor 
distributions  can  be  eliminated  from  these  transport  equation  results.  Going  from  the  far 
outer  gas  edge  toward  the  metal  surface,  whose  location  has  yet  not  been  determined,  the 
following  sequencing  logic  should  be  adopted.  For  each  nodal  point,  we  need  to  determine 
whether  liquid  metal  existed  there  at  the  previous  time  step.  It  is  assumed  that  once  the 
gas  phase  appears  it  will  remain  gas  and  liquid  metal  cannot  appear  again  at  that  nodal 
point.  If  gas  phase  existed,  x  will  be  remain  equal  to  unity,  and  the  temperature  will  be 
computed  from  the  calculated  local  enthalpy  and  gas  species  mass  fraction.  This  assump- 
tion is  justified  by  the  fact  that  in  most  combustion  stages,  metal  vapor  is  totally  consumed 
in  the  reaction  zone,  and  its  concentration  during  the  ignition  period  is  very  small. 

By  assuming  a  constant  value  of  C  ,  enthalpy  for  each  species  can  be  defined  as  fol- 
lows : 

For  fully  melted  metal  ingredient  at  T^  temperature. 
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hMl    =    CP„  (  W  298-  '    +    hH,f  ^ 

Ms 


For  fully  vaporized  metal  ingredient  at  "T-. 


Mb- 


\e  ~  h«i  +  cPMf<  TMb  -  T«,  >  +  hHfg  <6-10> 


Enthalpy  of  vapor  metal  ingredient  when  T  is  greater  than  T_   is  given  by: 


V,   =    hM3    +    Sj  T  -  TM..  >  <6"11' 


Enthalpy  of  steam  ingredient  at  the  saturation  temperature  will  include  its  heat  of  for- 

0 

mation,  h^  ,  at  298.  K  and  phase  change  latent  heat,  hwf  . 


h        =  h  °  +  C      (T  .  -298.  )  +  h  .  (6-12) 

wsv  w  p    ,v      wb  '  wfg  * 


Water  vapor  ingredient  enthalpy  at  temperature  above  the  saturation  point  will  include 


hwsv 
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h      =  h        +c      (T-T.)  (6-13) 

wv  wsv  p      x  wb  '  \  ' 

wv 


Hydrogen  product  of  combustion  will  always  appear  in  the  gas  phase,  and  its  enthalpy 
term  will  be  : 


hH  =  cp  (  T  -   298.  )  (6-14) 


PH 


According  to  Eq.  6-2  the  total  enthalpy  of  the  gaseous  mixture  is  defined  as  : 


h  =  hMcM  +  h    c      +  hucu  (6-15) 

M    M  wv   wv  H    H  v  ' 


After  introducing  the  terms  for  each  ingredient  enthalpy  Eqs.  6-11,13,14  and  defining  C 

g 
as  : 


C      =cMcM  +  c      c      +ccu  (6-16) 

Pg  pMg   M  P^wv  ph    H 


The  temperature  at  the  nodal  gas  phase  can  be  eliminated,  and  its  expression  is  : 


T  -  th  -  <  h«3  -  cp„  TMb  >c«  -  <hw-  -  cp    373->c^  "  cphch298-1  '  (  C      )  (6-17) 

Ma  wv  n  g 


Mg 

22 


If   at  the  previous  time  step,  x   was  less  than  1,  or  even  equal  to  0.,  the  procedure  for 
determining  a  new  value  for  x,  should  be  as  follows  : 

By  knowing  the  local  metal  vapor  mass  fraction  and  the  total  local  molecular  weight, 
assume  the  local  temperature  to  be  equal  to  T^ 


T^  =  A/[B-£n((cMMt)/MM)]  (6-18) 


For  aluminum  A=  36560.8  and  B=  25.016. 

Metal  vapor  quality  value,  x,  can  be  eliminated  from  the  total  enthalpy  expression.  By 
assuming  that  the  nodal  ingredient  mixture  consists  of  saturated  metal  vapor,  water  vapor 
and  hydrogen,  and  assign  the  local  temperature  equal  to  T_  x  remains  the  only  unknown 
and,  therefore,  can  be  eliminated: 


x  =  [h  -  (hMl  +  c      (T„.  -  TM  )cM  -  (h        +  c      fT-  -  373.))c 

Ml  p„,v    Mb         Mm/    M      v    wsv  p      x    Mb  "    wv 

Mt  wv 

"  CpJTMb  -  298>c J'  <  hMfgCM)  <6"19> 


If  the  resulting  value  of  x,  based  on  those  assumptions,  is  larger  than  one,  then  the  local 
temperature  should  be  greater  than  the  assumed  temperature  "T- ,  and  a  new  tentative 
value  of  temperature  should  be  determined  by  using  Eq.  6-17.  If  ,  at  the  first  iteration  x 
becomes  equal  to  or  greater  then  zero,  then  the  temperature  is  set  equal  to  T^  and  the 
surface  location  is  determined.  Further  inside  the  metal,  temperatures  are  determining  by 
setting  x  equal  to  zero. 


T=(h"hMl   +cp    TmJ'   cp„,  (6-20> 


PMf   ""  PMf 
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The  thermodynamic  relationship  for  the  metal  vapor  density  shows  that  increasing  of  x  will 
decrease  the  density  pM,  Eq.  6-21.  Metal  surface  vaporization  cause  a  local  density 
decrease  and  since  the  continuity  law  must  also  be  conserved,  the  net  metal  vapor  will 
expand  outward  to  the  gas  surrounding. 


P„=  PMg/[1   +  d  -  x  )(  pMg/pMf  -  1  )]  (6-21 ; 


6.2  Determining  the  Effective  Radiation  Heat  Source 

Considerable  amounts  of  radiating  heat  fluxes  govern  the  metal  combustion  process  in 
the  case  of  large  volatile  particles.  In  most  other  cases  the  surrounding  gas  can  be  con- 
sidered as  an  optically  thin  medium.  When  oxide  particles  appear,  the  media  is  not  trans- 
parent any  more.  The  number  of  particles  per  unit  volume  N,  the  particle  size  r  ,  specific 
particle  emissivity,  cloud  thickness  I,  and  the  relative  geometric  location  of  the  particle 
cloud  determines  the  radiation  fluxes. 

Intensity  attenuation  of  radiation  propagating  through  a  particle  cloud  is  given  in  Ref. 
13  as: 


l/lo  =  EXP(  -el  )  (6-22) 


Where  the  extinction  coefficient  e  is  a  function  of  N  and  r  .    Measurement  results  described 

p 

in  Ref.  13  demonstrate  the  validitiy  of  the  relationship  between  these  parameters.  We 
present  in  Fig.  5  some  data  results  of  aluminum  oxide  cloud,  which  will  be  used  as  a 
data  base  for  our  model. 
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Fig.  5  Correlation  of  all  spectal  data.  ( from  Ref.  13  ) 
Solution  of  the  gas  phase  chemical  reaction  along  with  the  other  conservation  equa- 
tions will  produce  the  oxide  product  distribution  in  the  field.    By  assuming  an  average  par- 
ticle size  r  ,  and  knowing  the  local  cell  density,  we  are  able  to  determine  the  number  of 
local  oxide  particles  in  that  particular  volume  cell  : 


N.   m  (3./4.KI-    )  (p./pj  =  <%2V/(4rcr  &S"3./3) 


(6-23) 
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By  accounting  for  ail  nodal  N.'s  and  their  radial  position,  it  is  simple  to  find  the  total  parti- 
cles cloud  depth  and  its  average  number  of  particle  per  unit  volume. 


N  =  (  LN.AV.  )  /  (  LAV.)  (6-24) 


Using  the  chart  in  Fig.  5,  we  can  find  the  value  of  e  then  the  overall  cloud  emissivity  e  ,  is 
determined  from  : 

zp  =   1  -  EXP  (  -  el  )  (6-25) 


Fig.  6  shows  directions  and  radiation  fluxes  interchanges  between  the  cloud,  the  metal 
surface  and  the  environmental  gas.  For  practical  purposes  it  is  assumed  that  all  media 
involved  are  grey.  As  the  gas  is  optically  thin,  and  particles  are  mainly  concentrated  in  a 
narrow  zone  (usually  the  flame  zone),  isothermal  radiating  surfaces  is  a  resonable 
assumption.  Total  particle  cloud  radiation  is  diffused  toward  the  metal  surface  and  to  the 
outer  gas  space.  Part  of  the  radiation  emitted  from  the  cloud  and  impiging  on  the  metal 
surface,  will  be  reflected  back  to  the  cloud  by  the  grey  metal  wall  and  absorbed  either  by 
the  particles  in  the  cioud  or  by  the  gas.  The  net  result  of  radiation  fluxes  through  the  par- 
ticles cloud  is  also  a  function  of  the  mean  gas  and  particle  cloud  temperatures,  T  ,  T  and 
wall  metal  temperature  T 


(q  "/ff)  =  atT4F      s.   +  a  c  T  %     +  a  (e  T    4(1-a  )F      S.-eT    *         (6-26) 

x^ox  '         p  w    w      w-c    i  p  g    g      o  pv  p    ox  x        w'    w-c    1       p    ox  v 

net 


The  s.  and  s    are  parameters  of  surface  ratios,  as  shown  in  Fig.  7  and  for  spherical  geom- 
etry they  are  given  by  : 


s.  =  A  .  /  (  A  .  +  A      )  =  r./(r.     +  r       )  (6-27) 

1  CI         V         CI  CO   '  CI  V    CI  CO      '  y  ' 
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s=A/(A.+A     )  =  r      /(r  .     +  r       )  (6-28) 

O  CO         v         CI  CO   '  CO       y    CI  CO      '  v  ' 


The  geometrical  view  factor    F       is  given  as: 

»  c-w        a 


F        =  [1./  (e  "'  +    (e  -1     -  )(  A  /  A  .))]  (6-29) 

c-w  L  v   w  v    p       -   1    /v       W         Cl//J  v  ' 


The  metal  surface  absorbs  radiation  from  the  particle  cloud  and  the  gas  while  emitting  its 
radiation  heat  flux  outward. 


(q       *i~\  -  „  *  t    4F        +oeT4-eT4  (6-30) 

v^w      .   /0)  —    a   E     I  c-w  w  g    g  w    w  x  ' 

net       '  w  p    ox 


Those  two  radiating  fluxes  divided  by  the  appropriate  characteristic  length,  which  is  the 
particle  cloud  thickness,  provide  the  local  corresponding  heat  source  that  should  be  taken 
into  account  as  source  terms  in  the  energy  equation. 
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Fig.  6  Radiation  heat  fluxes  exchange  between  the  oxide  cloud 
and  the  metal  surface. 
6.3  Energy  Equation  for  the  Oxide  Particle  Phase 

The  smoke  oxide  product  of  a  volatile  metal  combustion  requires  an  energy  equation, 
as  it  is  a  second  phase  in  the  field.  Since  the  small  oxide  particles  are  dispersed  in  the 
gas,  conduction  heat  exchange  between  them  can  be  neglected.  The  heat  transfer  time 
constant  for  the  largest  oxide  particle  (which  is  less  than  2  micron),  in  a  laminar  flow  is 
less  than  10  microseconds,  therefore  the  radiation  heat  losses  of  the  particle  cloud  may  be 
incorporated  as  a  source  term  in  the  gas  phase  energy  equation.  If  local  temperature  is 
smaller  than  the  OBT,  then  immediate  condensation  of  all  the  new  amount  of  the  oxide 
product  will  occur. 
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Fig.  7  Geometrical  relationship  between  the  cloud  and 
metal  droplet. 
If  M2  represents  the  volume  fraction  of  phase  2,  then  the  contact  surface  A     between 
the  oxide  particles  with  an  average  radius  r    in  a  volume  V  is  given  by  Ref.  14  as: 


A    =  3^2V  /  r 

s  P 


(6-31) 
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The  source  term  for  that  phase  will  include  internal-phase  heat  transfer  and  release  of 
latent  heat.  Ref.  14  presents  the  source  term  as  : 


Q2     =  Xc(As/V)(  Ta  -  Ti  )  -  Wqx  hoxfg  (6-32) 


The  energy  equation  will,  therefore,  be  presented  as 


d(p2h2)/d\.  =  Q2   -  p2u2dh2/dr  (6-33) 


The  oxide  temperature  can  be  directly  found  from  its  enthalpy  h2 


Tox  =  (h2-hox°-Cp    298-)  '  CP  (6"33a) 

ox  ox 


6.4  Mass  Transfer  Equation  for  Each  Species 

Analytical  descriptions  of  species  mass  fraction  and  enthalpy  transport  are  of  the  same 
kind  as  the  general  diffusion  equation.  At  some  time,  when  Lewis  Number  equals  1.,  their 
non-dimensional  solutions  are  the  same.  The  distribution  in  time  and  space  of  both  vari- 
ables, c.  and  h,  are  given  by  the  solution  of  a  diffusion  equation.  With  the  assumptions  of 
quasi-steady  combustion  and  unit  Lewis  Number  unity,  the  numerical  problem  reduces  to 
an  algebraic  set  of  9  equations,  (Ref.  8).  But  in  many  particular  combustion  cases,  such  as 
an  exploding  wire  or  a  composite  propellant,  high  temperatures  cause  dissociation,  so  that 
assumption  is  not  valid  any  more.  Simultaneous  solution  of  the  diffusion  equation  is 
needed  with  careful  attention  to  various  source  terms  and  temperature  dependence  of  the 
species  properties.  The  source  term  in  the  mass  transfer  equation  is  the  mass  production 
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or  consumption  of  the  each  particular  species.  At  each  time  step  and  location,  the  source 
term  is  determined  by  local  thermodynamic  condition.  If  the  local  temperature  is  not  high, 
enough  the  species  source  term  will  be  determined  by  the  reaction  rate   as: 


i)        i) 

WM"  _    R  TnFXPf-F  /R-nCM  mCu«  u^  (6-34) 


'm    =  B  -TEXPl-E^RTrM  mUwv  wv 


If  the  reaction  rate  consumes  more  than  the  species  mass  flux  into  the  local  volume  cell, 
then  the  local  consumption  and  production  are  controlled  by  the  'fuel'  and  'oxidizer'  diffu- 
sion rates.  If  the  local  mixture  ratio  f,  Eg.  6-35a  is  less  than  the  stochiometric  mixture  ratio 
fS{  (Eq.  35b),  then  the  new  value  of  the  local  oxidizer,  (water  vapor),  mass  fraction,  will  be 
governed  by  the  mixture  fraction  and  all  the  fuel,  (metal  vapor),  will  be  consumed.  A  new 
value  for  fuel  mass  fraction  will  be  set  and  complete  consumption  of  the  oxidizer  will  occur 
when  f  is  greater  than  f&±.  This  procedure  is  presented  graphically  in  Fig.  8,  and  Eqs.  36 
show  its  analytical  representation. 

f  -  1./(1+cwtf/cM)  (6-35a) 

f8t  -  1./(1.  +  (cJcH)    )  (6-35b) 

st 


For   0<f%ofst  c  =  0  (6-36a) 

M 

C       =    1.-f/f    . 

WV  ST 

For      1      f$fst  cm  =  0.  (6-36b) 

cM  =  (f-C)/(1-fst) 


The  products  mass  fractions  are  determined  according  to  the  stochiometric  ratios  i.  One 
Kgm  of  aluminum  will  produce  iQx=  1.88653  kgm  of  alumina  and  iH  =  0.1 1898  kgm  of  hydro- 
gen, in  the  process  it  will  consume  i     =  0.998  kgm  of  water.    The  net  rate  of  oxide  produc- 
tion will  depend  on  the  mixture  condition  and  mass  fraction  of  the  reactants. 
For  a  fuel  rich  mixture  the  oxide  mass  fraction  source  is: 
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Fig.  8  Mass  fraction  distribution  with  mixture  fraction. 


w       =  (i     /  i    )d(cMpi^fci)/dt 

ox  x  ox       wv'    x    It  ' 


(6-37a) 


and,  for  the  fuel  lean  case,  it  is: 


ox  ox    v    Mr  ' 


(6-37b) 
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Since  the  hydrogen  product  will  always  appear  in  the  gas  phase,  we  are  allowed  to  set  its 
mass  fraction  equal  to: 


cu  =  1-c  -c  (6-38) 

H  M     wv  v  ' 


The  diffusion  equation  for  the  mass  transport  of  species  j  is  given  by; 


.2         2 


d(pxc.)/dt  =  r   d[r  PlD.o{dc./dr)]/dr  +  w.    -  p^^dc./dr)  (6-39) 


The  density  is  the  density  of  the  local  mixture,  and  Di0  is  the  species  diffusion  to  the 
mixture.  Late  in  chapter  6.7,  we  will  show  how  to  specify  the  local  properties. 

6.5  The  momentum  equations  for  each  phase 

The  two  phases  have  the  same  type  of  the  momentum  equation  with  a  source  term 
that  stands  for  the  momentum  exchange  between  them.  We  assume  the  second  phase 
consists  of  solely  oxide  particles.  Since  the  particles  are  dispersed  in  the  gas,  they  do  not 
exert  friction  forces  on  each  other,  and,  therefore,  we  may  ignore  the  momentum  transport 
within  that  phase.  We  now  define  the  source  term  of  momentum  interphase  transfer.  From 
Stock's  flow  around  one  particle  we  know  that  the  friction  force, F    ,is  given  by  : 


Fp  =  6r^rp(  u,  -  u2)  (6-40) 


Now  we  define  F      as  the  force  F    per  unit 
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particle  volume. 


F*-  F/(4nr3/3)  =    4.5  fi/r*  (6-41) 


If  0t\  and  Ski  are  the  volume  fractions  of  phase  1  and  2  respectively,  then  the  momentum 
interphase  exchange  per  unit  volume  and  time  will  be  presented  as  follows: 


M12  =  -  M21   =  Fp  9Lx0t2  V(  Ui  -  u2  )  (6-42) 


The  radial  one-dimensional  momentum  equation  for  each  phase  will  then  be  given  by: 


pi(0Ui/ft)  =  r"2d  [r2p,n,(<?u,/dr)]/dr  -  p,u,(aii/dr)  +  M21  (6-43) 

p2{du2/dt)  =  -  p2U2(5u2/^r)    +  M12  (6-44) 


Since  we  solve  the  gas  and  the  metal  as  one  phase  by  ignoring  surface  tension  forces 
between  them,  the  metal  droplet  may  be  smeared  into  the  gas.  In  order  to  prevent  that, 
both  phase  velocities  inside  the  metal  should  be  set  equal  to  zero.  For  small  enough  time 
step  the  procedure  of  freezing  velocities  will  capture  the  oxide  on  the  metal  surface,  and 
therefore  simulate  deposition. 
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The  total  diffusion  of  species  j  into  a  mixture  which  contains  a  mole  fraction  x  of  that  spec- 
ies, is  defined  by  D.Q: 


D j0  =  (  1  -  Xj  )  /  (S(Xk  /  D  .R))  (6-49) 


The  mole  fraction  %.,  can  be  found  from  c.  ,  M^  and  M.  as  follows: 


x.  =CjMt/M.  (6-50) 


The  mixture  conductivity  and  viscosity  coefficients  can  be  found  from  gas  tables  Ref.15. 
They  are  expressed  as  a  fourth  degree  temperature  polynomial  for  each  species  present  in 
the  mixture.  The  conductivity  k.  ,  the  viscosity  \i.  and  heat  capacity  c      will  be  presented 


as: 


k.-Ek-.T1  (6"51a) 

H,-  ^T1  (6"51b) 


c     =  Ec      T 


l  (6-51  c) 
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In  order  to  evaluate  the  mixture  viscosity,  a  bi-molecular  viscosity  parameter  £1..  should  be 
defined  : 


ft3k=(1./V8)(1   +  M./Mk)'l/2[1   +  (ix.  I  jik),/2(  Mk/M.  )V4  ]2  (6-52) 


H  is  defined  as: 


H  =  I  [x3H3  /  (  2(xk^jk  )  )  ]  (6-53) 


The  mixture  conductivity  will  be  written  as: 


K-  I  [Xjkj  /  (  S(xk"3k  )  )  ]  (6-54) 


The  density  pi  is  defined  by  the  ideal  gas  equation  of  state. 


Pi  =  p/(RTi)  (6-55) 


c       was  defined  previously,  in  6-1  when  discussed  the  gas  energy  equation.   The  next  step 

PMg 

should  be  to  cast  the  model,  presented  above,  in  the  form  of  a  computer  program. 
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7.  SUMMARY 

Preliminary  calculation  for  aluminum  particles  from  a  shaped  charge  revealed  that  the 
maximum  size  of  a  particle  that  may  reach  ignition  conditions  while  moving  at  speed  of 
7000.  m/s  during  200  microseconds  in  a  steam  environment,  is  just  6  micron. 

Exploding  wire  experiments  in  water  vapor  chamber  is  an  adequate  technique  for  find- 
ing the  most  combustion  effective  metal  liner  of  underwater  shaped  charge.  However,  it 
does  not  simulate  the  phenomena  of  ignition  and  combustion  of  high  speed  moving  metal 
particles  in  the  U/W  S/C.  The  need  for  precise  prediction  in  any  of  these  cases  requires  a 
solution,  which  will  account  for  the  transient  effects  along  with  two  phase  flow  and  radia- 
tion fluxes  emitted  from  the  oxide  cloud.  The  combustion  process  parameters  that  will 
generalize  the  model  for  both  the  volatile  and  non-volatile  metals  cases  are:  the  metal 
vapor  quality,  surface  temperature  and  metal  vapor  surface  concentration.  By  writing  the 
energy  equation  for  total  enthalpy  that  include  those  parameters  and  accounting  for  local 
metal  and  oxide  phase  change  temperatures,  the  solution  will  be  valid  for  both  type  of 
metals.  The  conservation  equations  were  written  for  the  entire  region,  gas  and  liquid.  The 
technique  of  iterating  in  order  to  find  the  metal  surface  was  described.  Radiation  is 
accounted  for  as  source  terms  in  the  oxide  and  at  the  metal  surface. 

Such  a  solution  for  a  general  type  of  metal  burning  will  simplify  the  computer  code 
program  for  the  various  cases. 
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